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Heavy-to-heavy quark decays at NNLO 

Alexey Pak and Andrzej Czarnecki 
Department of Physics, University of Alberta, Edmonton, AB, T6G 2G7, Canada 

Details of a recent calculation of ©(t^s) corrections to the decay b — > c£vi, taking into account 
the c-quark mass, are described. Construction of the expansion in the mass ratio mc/m;, as well as 
the evaluation of new four-loop master integrals are presented. The same techniques are applicable 
to the muon decay, ^ eVyM^. Analytical results are presented, for the physical cases as well as for 
a model with purely-vector couplings. 



INTRODUCTION 



Because of the role that the semileptonic decay 6 clvi and the muon decay /i 



play in the determination 



of the parameters of the Standard Model, it is warranted to improve the theoretical description of their rates and 
distributions. Recently, we have determined 0{^o?^ corrections (a^ for the muon), including the effects of the charm 
quark (electron) mass pj] . Since this was the first analytical evaluation of such mass effects and required an extension 
of known results, in this paper we describe some technical details. 

We made extensive use of earlier results obtained in the massless case 0, B [3| , in particular of the master integrals 
determined in those projects. It turned out, however, that the massive case requires the knowledge of more terms 
of the expansion of those integrals in the dimensional regularization parameter e — ^-j^. To find them, we took a 
different approach to the evaluation of those master integrals. That approach is described below, together with the list 
of new terms that are now known. Since the decays considered here are a model for other beta-decay-like processes, 
those integrals will Hkely find other appHcations in the future. 

Before presenting that calculation, in the following Section we introduce the notation and in Section IIIII describe 
how the expansion around the massless case is constructed. Our results for the corrections to the 6 — > c decay and to 
a decay in a model with pure vector coupHngs, as well as for the new terms in master integrals, are collected in three 
appendices. 



II. DECAY RATE AND RADIATIVE CORRECTIONS 



The tree-level decay rate and the one-loop corrections are known exactly Q- Their expansion in p = mc/rrib ^ 1 
and parameterization of two loop corrections are 
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Ca = 3. Here Nl ^ 3 



is the number of massless quarks, and Nh,c — 1 label, respectively, the 6-quark loop, and virtual and additional real 
c-quark contributions; components Xl, Xh, Xq, Xa-, and X]sia are separately gauge-invariant and finite. The limit of 
these functions for mc = is known [3!|. The purpose of this paper is to obtain their mass dependence. All calculations 
are performed in D = 4 — 2e dimensions, and axial currents are treated according to the prescription of . General 
gluon gauge was used to ensure gauge invariance of the results. 

The presence of the logarithm of the mass ratio in the lowest-order rate, Eq. ([2]), signals the presence of higher 
powers of logarithms in the higher-order terms (indeed, we see a quadratic logarithm in Eq. ([3])). It is for this 
reason that the term p"* requires more terms of the expansion of Feynman integrals in powers of e mentioned in the 
Introduction. 

Using the optical theorem, we find X2 as the imaginary part of 39 self-energy diagrams such as the examples shown 
in Fig. [TJ Each diagram is expanded in asymptotic regions 0, S, [§| to several orders in p and all contributions are 
summed. Some details of that expansion are described in the next section. 
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Figure 1: Examples of ©(a^) b-quark self-energy diagrams whose cuts describe the semi-leptonic decay. Dashed lines represent 
a charged lepton and a neutrino, whose masses we neglect. 




Figure 2: Expansion of a double-scale topology (a) in all contributing asymptotic regions (b-1). Thick lines represent mass nib, 
thin - mass rric, dashed lines are massless, double lines correspond to eikonal (static) propagators. Different regions correspond 
to expansions with different assignment of scales ("hard", or "soft", nic) to the four loop momenta. 



III. ASYMPTOTIC EXPANSION 



As an example, we consider the rightmost diagram in Fig. [T] which has the richest structure of asymptotic regions. 
Evaluating the traces and Wick rotating, we obtain a number of terms with the following factors in denominators: 

Di = e, D2^{q- r)2, Di^q^ + ml, D4 = (q + if + ml (5) 
D5 = {q + l-rf + ml , = = ^^ Ds ^ + 2pl, Dg = (p + k - q + rf, 

where fc, q, r, and I are loop momenta, and the external momentum p is on shell (p^ = —ml) (Fig. [2l[a)). Each loop 
momentum is then assigned one of the two relevant scales - mb ("hard") or mc ("soft") , and integrands are expanded 
in each case. All contributing regions for this topology are sketched in Fig. (Hb-l). As an example, consider q and I 
being "soft", and k and r "hard". The denominators are Taylor expanded as 



' and similarly (6) 



1 _ 1 
Ih ^ /2 _^ 2pl 2pl ^ V 

[(p + fc + r)2r' 

producing, together with the remaining denominators, the case (i). 

As one can see, in most regions this topology factorizes into known one- and two-loop integrals. Note the unusual 
"eikonal" denominator l/{2pl + iO) in the above example and similar factors in case (1), which lead to odd powers 
of p in the end result. Disentangling of the products of loop momenta in the numerator, naturally appearing in the 
expansion, in the most difficult cases involves the algebraic method of Ref. jlB| and, at high orders in p, leads to a 
huge number of terms. With intermediate expression size reaching hundreds of gigabytes, this calculation was possible 
only with the computer algebra system jlll] . 
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The three- loop eikonal i nteg rals in region (1) correspond to topologies studied in Ref. Using the integration- 
by-parts identities [H, we reduce the expressions to a few "master integrals" found in Refs.fll. and a 
previously unpubHshed integral calculated by V. A. Smirnov, Eq. IjCip . 

The four- loop "all-hard" case (b) here is the most challenging. To expand X2 to 0(^p^), an implementation of 
algorithm [I^ was running for several weeks. Finally, all diagrams were reduced to the 33 master integrals. Their 
evaluation was the biggest challenge of this work. 



IV. EVALUATION OF FOUR-LOOP MASTER INTEGRALS 

The results for master integrals given in [J| are sufficient to obtain 0[p^) and 0[p'^) terms of X2. However, in 
order to find the following 0[p'^) contribution, many integrals need to be expanded further. In j^l, the initial terms of 
that expansion were obtained in a series of steps. First, some of the internal lines of the diagram, representing a given 
master integral, were assigned a mass M much larger than the mass of the external particle m. Next, the diagram 
was expanded in m/M using a similar approach as described in the previous section; several terms of that expansion 
were obtained in each order in e. Finally, and most non-trivially, a pattern in that expansion was recognized and the 
expansion was now summed analytically. The value of the result at M = to gave the desired value of the integral. 

In higher orders in e, the recognition of the expansion pattern may be very difficult. Instead, we evaluate these 
integrals using the method of differential equations We start by choosing a few artificial double-scale topologies, 
related in some Hmit to th e ne eded integrals. We illustrate this approach with the topology shown in Fig. [31 It 
involves an artificial mass y^l/x, and in the on-shell limit a; ^ 1 reproduces the topology needed for the first diagram 

shown in Fig. [1] (where we also introduce the loop factor = ^^^2)- 

This function is chosen due to several reasons. First, if < a; < 1, the cuts are the same as in the on-shell limit, and 
we may discard the real part in all calculations. Second, large mass expansion to a few orders in x and e is relatively 
simple Q- And third, the associated differential equations have a structure convenient for an iterative solution; this 
property will be explained later. 

This topology has 40 master integrals. Derivative of any such integral can be re-cast in terms of integrals with 
shifted indices. 



da d 

— /(ai, ...,ag;x) = ■■■'°^5 + 1, •■•) + :;2'^(°^i' ■•■'°^6 + 1, ■■■), (7) 



generating a system of 40 differential equations. It can be split into independent subsystems of at most four equations, 
where the solution of one system enters the RHS of the following one. 

For example, consider functions /(x, e) = /(0, 1, 1, 0, 1, 0, 0, 1, 0; x), g{x,e) = /(— 1, 1, 1, 0, 1, 0, 0, 1, 0; x), and 
h{x, e) — 1(0, 1, 1, —1, 1, 0, 0, 1, 0; x), entering a particularly simple system of relations 

4x- 5 + 6(7 -6x) , , 9 -lie 3e - 3 ^ 

f = 2x{l X) ^+2(1^-'^ +2(1^^' 

1-e 8 + x-e{10 + x) ^ - 3e 

^ ~ 2x{l-x)-^ 2x(l - x) 2(1 -x) ' 

,/ 1 — 1 — e 3e — 3, 

We solve this system with respect to derivatives of /, and expand f{x,e) — jf-i{x) + fo{x) + efi{x) + .... Finally, 



1 f [d^k][d^q\[d^r] D-^^ 
/(ai, ag; a;) = — -=^lm I ^ , 

^ ^ 7rJF4 j D^' D^'' D^<' D^-' D^^ 

= -l,I>i = k^, D2 = {k + pf, D3 = {k + q + pf, 
D4 = {k + q + r + pf, D5 = {q + r + pf + 1/x, 
Da = {r+pf + 1/x, D7 = g^ Ds = r^ Dg = 2qp. 



Figure 3: Auxiliary double-scale topology. Thick lines represent mass ^JX/x, thin - unit mass propagators, and e labels a 
denominator raised to a non-integer power. 
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the differential equations for fi become 

„ 9-6x 18 -6x 6 

^ x{l-x)^' x^l-x)^'^ x^il-x)-''"^'' 
fi. -n fi- - IQ - 47- 9x , 31 

The three solutions of the homogeneous equation (with Ri = 0) can be guessed: 1/x, 1/x'^, and (1 — x'^{x + 6) + 
6x{l + x) lna;)/a;3. Euler's formula allows then to solve the inhomogeneous equations. To fix the three integration 
constants, we use the large mass expansion 

, 1 1 15 + 4 In a; , , 

= -4^+24 8^ ^ ^ 

2Tr2 — 145 — 601na; — 81n a; 1 In a; x x"^ 



16a; 4 12 144 1440 

The general solution for fi is expanded and matched to this series to find suitable constants. Finally, we obtain 
a;-dependent solutions: 

, _ 1 ,1 15 + 41na: .19 1 67r2 - 373 

- 'Ji' ^° " 24 8^ ' " 36 + 12^ + 48.x ^^^^ 

ff(0;x)(x-45) ff(l;x)(x3 + 9x2 -9x- 1) iy(o, 1; x)(l + x) iJ(0,0;x) 

12^^^ 12^3 2x2 X ' ■ 

Due to the fact that the answers can be expressed in terms of harmonic polylogarithms (HPLs) H(...;x) fis'l, the 
expansion can be continued as long as CPU resources allow. (More accurately, this procedure works as long as the 
solutions of the homogeneous equation, as well as their inverse Wronskian and minors of Wronski matrix depend on 
X in the denominators only through factors x and 1 ±x, and numerators contain HPLs and polynomials. This allows 
to solve the inhomogeneous equation in terms of HPLs.) Taking the limit x — > 1, we obtain the required on-shell 
integral. 

One integral of the topology shown in Fig. [31 m(x, e) — /(0, 1, 1, 1, 1, 1, 1, 1, 0; x), presents an additional difficulty. 
Its x-dependent expression starts at 0{e) and logarithmically diverges at x = 1. To find that integral, we evaluated 
instead the infrared-safe function /(1, 1, 1, 1, 1, 1, 1, 1, —1; x), took the on-shell limit and then reduced to m(1, e). This 

lead to finite 0(e°) and 0{e) expressions: u(l,e) = - if^ - e + ^^f^ + ff^)- As a tradeoff, other integrals 
were needed to a higher order in e, involving harmo nic p olylogarithms up to weight six. 

An important tool was the software package HPL and additional software for the series expansion of harmonic 
polylogarithms was developed [20]. As an independent check, we used numerical integration of the Mellin-Barnes 
representation of some integrals [21]. 



V. RESULTS AND SUMMARY 

The rather lengthy expressions for Xl, Xh, Xc, Xa, and Xna calculated through 0[p^) are presented in Appendix 
[Al Their features, such as the logarithms and even and odd powers of p, together with physics consequences, have 
been discussed in Here we would like to illustrate the convergence of the expansion. To this end. Fig. |4] shows 
the plots of Xi as functions of p, in two versions: soHd lines show all known terms while dashed fines are obtained by 
leaving out the last known power of p. We see that the convergence is excellent up to at least p — 0.3, of interest in 
this study. 

In addition to the integrated decay rate, precision fits to experimental data are done for the moments of the lepton 
energy Ei and the hadronic-system energy Eh distributions in the rest frame of the 5-quark, with the goal of accurately 
measuring several parameters including iVcbl, mf,, and Wilson coefficient of non-perturbative operators. Thus, QCD 
corrections to those moments are also of interest. These corrections are defined by 
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and similarly for the moments of Eh , described by coefficients Hj ; the average is taken over the whole phase space 
of decay products. One-loop spectra are known exactly [l^l, and Fig. [5] presents NNLO corrections L^2'^^ ^^'^ E[^''^\ 
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Figure 4: Mass-dependent corrections to X2 of Eq. ([TJ. 




Figure 5: First two moments of lepton and hadron energy Figure 6: Charm quark contribution to semi-leptonic 6 

distributions u decays 

In [H we discussed how the experimental cuts can be approximately modeled by the analytical calculation. Finally, 
we present the analytical results for Uc (Eq. I|A6P and Fig. [6]), an analogue of Xq in Eq. ([T|), describing charm quark 
loop contribution to the process with a massless quark in the final state, b — > utv. 

The same method can be used in the model in which chiral weak coupling of quarks, |^7/i (1 — 75), are replaced 

with a pure vector vertex, ^^7^^. This is a useful toy model, e.g., for logarithmic re-summation studies [23|- In 
parameterization similar to Eq. ([T]), the second order components Vi are defined in Eqs. (|Bll4B8p . We evaluated terms 
to and it is now easy to find more. 

To summarize, in the process of this calculation, we checked and confirmed the massless limit of the 0{al) cor- 
rections 0, 0; We extended those results to several orders in the mass ratio of the final and initial quarks. To 
complete this task, additional terms of master integrals were required, and the differential-equation method was ap- 
plied to compute them. Our results agree excellently with the numerical calculation jH]. An ultimate test of the 
convergence of the mass expansion will be a corresponding expansion around the opposite mass limit, rric — mt. Work 
on this is in progress. 
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Appendix A: ANALYTICAL RESULTS 

Here we present the components of the total decay rate evaluated to O (/?''). In the case of the 6-quark loop 
contribution, X^^ Fig. [4] shows an extremum point around p = 0.2. This contrasts with the other contributions, 
that seem to be monotonous functions of p. In order to double-check the convergence, also the 0{p^) term has been 
computed. The final results are 
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Charm quark contribution to 6 ^ uiv decays through 0(^p^) is: 
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Appendix B: A VECTOR MODEL 

As a byproduct of this project, for the purpose of comparisons with [l^l, we have determined 0(a1^ corrections to 
the decay rate in a model where the VF-boson has only a vector, and no axial-vector, coupling to fermions. That is, 
its interaction vertex is obtained from the standard {V — A) one by the substitution ^^7^ (1 — 75) ^7m- The 
results for this model are parameterized in analogy with Eqs. I|ll2l3l4p . with replacements Xt ^ Vi. Results for the 
tree-level and first-order corrections, expanded through ©(p^), are 
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The difference with the standard decay is seen already at the tree level: odd powers of p are present in Vb, while in 
the V — A decay they appeared only in the first-order corrections. The reason for this is a contribution of a four-quark 
operator, as discussed in [J, US]- In the V — A decay, a similar operator contributes only through its 0{as) matrix 
element. As a result, the second-order corrections in the vector case are even more complicated than in the V — A 
decay. They read 
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Appendix C: MASTER INTEGRALS 

In this appendix, we collect new results for master integrals that were obtained in this project. 
For the three-loop integrals needed for the factorized diagrams of the type shown in Fig. [2] (1), the 0{e) term of 
the integral J3 of Ref. was needed. It was privately communicated to us by V. A. Smirnov, 

J3 _ 32 / 256 In 2 448 647r\ ^_ r(l + e) 

^2^3 3^ 3 3 3 y ^' (47r)^/2 ' 

For the un-factorized four- loop diagrams in Fig. [Hb), the master integrals have already been classified and largely 
evaluated in Ref. Q. Using the approach described in Section [TVj we obtained additional terms of their expansion 
in e. To save space, we present here only those additional terms, using the notation of Ref. in particular Eqs. 
(5.7-11) of that paper. Instead of repeating here the long expressions already known, we define Ajy, N = A..F, to 
denote the new terms. For any affected master integral, its improved value can be found as follows, 

Im [In (fli, ■ ■ • , flu)] = 7rs^ (N^)"^ [terms already known + An (ai, • • ■ , an)] , (C2) 

where, as defined in Ref. 01, -^e = (^J)^ ^""^TXtScT^ * denotes the square of the momentum flowing into the 
diagram. Its power j depends on the integral considered and, for each one, can be found in Ref. 01. Table H] shows 
the extra terms A^r. 

We use the following notation for transcendental constants. An — Li„i, n — 4. .6, and se — 0.9874414.... Here Li„ 
denotes polylogarithms [l^l, ^'hiX = X^i^i W ^^"^ '^6 = (00) is certain harmonic sum 
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Table I: Additional terms for master integrals of Ref. f^, as defined in Eq. IIC2 



Aa 
Aa 
Aa 
Aa 
Aa 
Aa 
Aa 
Aa 
Aa 
As 
As 
As 

As 
Ab 
As 
As 



Ac 
Ac 



Ac 
Ae 
Af 
Af 
Af 

Af 
Af 



1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0) 55C5 + e - 2^ - 26^2^3 - 51C| + 220(5 

1. 1. 1. 0. 1, 0, 1, 1, 1, 0, 0) 6^ - iifP + ^ + ^ - i2z^ + 279C5) 

1 n 1 1 n 1 n 1 n 1 ITI e'^ I 2515231 _ 10153l7r^ _ 13097r^ _ 3595IC3 _ Ttt^Cs „ 779C5 \ 

±, U, 1, 1, U, ±, U, 1, U, ±, t I -^g^^ 2592 180 72 3 J 

1. 1. 1. 1, 0, 1, 0, 1, 0, 1, 0) 6^ lli^ + + 37^ + _ 41^^^ 

1 n 1 n 1 1 1 n 1 n n"! I 51175105 _ 2388197r^ _ 31497r^ _ 54895C3 i 437r^C3 _ 1243C5 

±, U, 1, U, ±, ±, 1, U, ±, U, t 1 15552 2592 360 72 3 2 

1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0) _6^ + 367^ + i2fC3_4^ + 59^^) 

1, 0, 1, 1, 1, 1, 0, 0, 1, 1,0) e^i^- llfp + ^+ 600C3 - + IO92C5) 

1, 0, 1, 1, 0, 1, 1, 1, 1, 0, 0) _ 147441 _ 14^ ^ 112P ^ 19,^ _^ ^gg^ ^ 

1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0) _ 38^ + 1^ + 2i«C^ _ 10i|!c^ + 563C5) 

1 1 1 1 1 1 1 1 1 1 n i245C5 „ ITzil _ TT^fCa 
J^i J^) J^) J^i J^i J^j J^j J^i J^i J^) '-J 32 360 32 

1,1,1,1,1,1,0,1,1,1,0) -1^^'-^^-^ 

1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0) 6^ (727r2 - 12f« + 2^ + 1^ - + SMCs) 

1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0) e (ifl - ^ - - 68C3 + Tv^s - 27(5) 

1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0) (2172 - Sllzd _ iSlgl „ 472C3 + - 232(5) 

1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0) 6 (24^4 + 6As~'-^-67vHn2 + ln^2-^ + 12z^ + 15(3 - + ^) 

1, 0, 1, 1, 1, 1, 0, 0, 2, 2, 0) 6^ (i^S _ 12^4 _ 6A5 - 4^ - _ ^2 ln2 - - i2z^ - + + ^) 

(12A4 - I2.I5 -6A6-i5|£a + l4l + i2ll2_l^_13|li + „ 17^2 1^ 2 - 



4 

975C3 j„ 477r^C3 _ 197r^ In 2 31l7r"^ In'' 2 5l7r''C3 In 2 , 721377^ , 2177C| _ 17490 



16 2 ' 3 40 160 ' 2 ' 5040 ' 16 



1,0,1,1,1,0,0,1,1,1,0) e^i£^»l<f5 _4^2i^2+i||f^ + 105C3) 
1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0) e (113 - 48^4 - 2 In'' 2 + + 4tt^ In 2 - + 14(3 



6 ' 60 

934 - 576^4 - 480^5 - 24 In'' 2 + 4 In^ 2 + + 427r^ In 2 - 87r^ In-" 2 + ■ 

r'* 3ir*ln2 , 24(3 — SHufil -|- Sl^Cs ^ 



72 2 



1,0,1,1,1,1,0,1,1,1,0) 6^ (i2i899 _ 927r2 In 2 + 192^4 + 8 In"* 2 + + IGn^ in" 2 - ^iil + ^ +7t^Cs- 45(5) 

1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0) 12z^ 94(5 + e (S22^ - 4|1 - 2^ + 6(1 - 376(5^ 



22l7r^C3 _ 799C5 
48 . 16 



1,1,1,1,1,1,1,1,1,1,1) ■ 

1,1,1,1,1,1,1,0,1,1,0) fiigCi - 3grl _ lOlSC^ 

1,0,1,1,1,1,1,0,1,2,0) ^e(^ + i^ + '-^ 



20 ' 48 

1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0) 6 (l031 - i2|Ei _ 22|e1 _ 163(3 + ^ - 64(5 j 

1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0) 6^ (25|!d _ 2172 + iOTni ggg^g _ 98^^ _|_ 4Q3^^^ 

+ ^3 ^2CW ^ 54^ _ ^7Q37 ^ 4Sg^ ^ 58gg^^ „ ggg^^l ^ ^^33^^ _ WTg^ ^ 

Af(1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0) 46yl4 - 30 + ^ilHll + 111 ^ 27v'Hn 2 - - ^ + 24(3 + e (184A4 - 272 + 46A5 + ■ 



3 



23 In^ 2 I 677r^ _ o 2 1„ 9 _ 73-;r'^ _ 207r^ In^ 2 _,_ 57r^ In^ 2 _ 5477r'^ In 2 , 209C3 _ 997r^C3 _i_ 3921C5 \ 
60 ~'~ 6 360 3 ~'~ 9 360 ^ 2 8 ~'~ 16 ^ 
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